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ABSTRACT 

A wind's turbine role is to convert the kinetic energy of wind into energy electric. Its components are designed to 
maximize the conversion. Different control structures have been developed according to the characteristic power speed. 
In our case, we propose the control strategy of the wind turbine to maximize power with control of the speed. 
This control strategy consists in determining the speed of the turbine which provides the maximum generated power. 
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INTRODUCTION 

During the last years, there was a strong penetration of the renewable resources of energy in the network of power 
supply. The wind power production played and will continue to play a very important role in this domain for years to 
come. 

Wind turbines have base of the doubly fed induction machine (DFIG) undoubtedly appeared as one of the high 
technologies for the manufacturers of wind turbines, Demonstrating that it is about an actual cost, Effective and a reliable 
solution. 

The first part is dedicated to the description and modelling of wind turbines based on physical equations 
responsive operation. 

The second part, we present a mathematical model of the (DFIG) will simulate the model in generator mode. 

The third is devoted to the study of the technique of indirect control power and compare the simulation results by 
PI Controller with Anticipation and PI Controller with Phase Advance. 

MODEL OF THE TURBINE 

The model is based on the characteristics of steady state power of the turbine (Lubosny. 2003). 

P m = j^Pmt = C p .P mt = ~C p m P nR 2 Vl (1) 

With A = — (2) 

v x 

n ± : Rotation speed before multiplier. 

R :rotor radius 35.25 m 

p: air density, 1.225 kg.rrT 3 
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C p = f(l fi) =C 1 ^f - C 3 p - C 4 ) exp g) + C 6 A 



(3) 



i i 
with : — — 



%i >,+0.O8P |r+i 



et C x = 0.5176; C 2 = 116; C 3 = 0.4;C 4 = 5;C 5 = 21; C 6 = 0.0068 [1]. 
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Figure 1: The Power Factor for Different Angles of Stalls 

Characteristics of C p in terms of X for different values of the pitch angle are shown in figure 1. 
The maximum value of C p (C pniax — 0.4353) is reached of P=2°and h= 10.01. This particular value of X is defined as the 
nominal value A nom (Lubosny. 2003; EL Aimani. 2004). 

DYNAMIC MODEL OF THE DOUBLY FED INDUCTION MACHINE 

A commonly used model for the doubly fed induction generator (DFIG) is the Park model. 
The electrical equations of the DFIG in the Park reference frame are given as follows (EL Aimani. 2004; Poitiers. 2003): 



iV sd = R s i sd + — -a) s (p. 



sq 



1 v sq = R s i sq + —jf + a> s (p sd 



n • i d<Prd 

I v rd — K r l rd + ~71 ^rtyrq 



'rq 



dqj rQ 



(4) 



(5) 



The stator and rotor flux are given as: 



(<Psd 
{<Psq 



Lshd + L m i rd 

sq L s i S q L m i r q 



(<Prd — L r i rd + Lmhd 
V^Prq L r l r q + L m i S q 



(6) 



(7) 



In these equations,R s , R r , L s and L r are respectively the resistances and the inductances of the stator and the rotor 
windings, L m is the mutual inductance. 

v sd< v sq< v rd > v rq< isd< 'sq' ird > hq> 'Psq' ^rd ' ^rq 316 ^ ^ anC ^ 1 com P onents of the Stator and rotor Voltages, 

currents and flux, whereas co r is the rotor speed in electrical degree. 



The electromagnetic torque is expressed as: 
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— PvPsd - hq ~ <Psq- hd) 



(8) 



Stator and rotor variables are both referred to the stator reference Park frame. With the following orientation, the d 
component of the stator flux is equal to the total flux whereas the q component of the stator flux is null figure. 2. 




-*^fJl \ 

Stat at axis 

Figure 2: Determination of the Electrical Angles in Park Reference Frame 

<Psd = <Ps, <Psq = 0 

By replacing (9) in (6) and (8), the electromagnetic torque can be given as follows: 
C 



(9) 



(10) 



Assuming that the resistance of the stator winding R s is neglected, and referring to the chosen reference frame, the 
voltage equations and the flux equations of the stator winding can be simplified in steady state as follows: 



f Vsd = 0 

hsq = V s = 0) s (p s 

(<Psd — L s i S d + L m i rd 
I 0 = L s i S q + L m l r q 

From (12), the equations linking the stator currents to the rotor currents are deduced below: 

r j — f s Lm ; 

"-sd ^ ^ *-rd 



(ID 

(12) 



(13) 



-sq 



The active and reactive powers at the stator side are defined as: 

(Ps = v sdhd + v sqi-sq 
IQs = v sqisd + v sdirq 

Taking into consideration the chosen reference frame, the above power equations can be written as follows: 

Qs = V s i S cl 

Replacing the stator currents by their expressions given in (15), the equations below are obtained: 
ft 



(14) 



(15) 



n _ Vs<Ps v s L m 
Vs — L ~ L l rd 



(16) 
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The block diagram of the DFIG model in Park reference frame is depicted in figure 3, assuming a constant stator 
voltage (vs). 
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Figure 3: Block Diagram of the DFIG Model 

REGULATION WITH BUCKLE OF POWER 

To improve the control system the DFIG, we will introduce an additional loop control of active and reactive power 
in the block diagram of the control loop without power so that each axis controller contains two PI control, one to control the 
power and the other rotor current (figure 4) (Mehdary. 2009). 
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Figure 4: Schema Block Indirect Regulation with Loop Power 
SYNTHESIS OF THE REGULATOR PI FOR THE CONTROL OF THE POWER 

It is a simple and easy to implement controller. Figure 5 shows a closed loop system corrects by a PI controller. 

k- 

In our case, the transfer function is in the form k p + — as shown in Figure 5 



R r + pL r c 



Figure 5: Diagram of the Control System hy PI 

The transfer function of the open loop (FTBO) with regulators is written as follows (EL Aimani. 2004; Moriarty, 
Butterfield. 2009; Multon, Obin, Gergaud, Ben Ahmed. 2003): 
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k i 



l 



+ V 



Applying the method of compensation pole transfer function 1, we obtain the following equality: 

kj _ gr 
k p L r a 

After compensation, the following are obtained FBO: 

bB0 ~ Way 

This gives us the transfer function following closed loop: 
F ~ — 

_ 1+T P 

With: t = ^ and a = (l - 

v. is the response time of the system that we fix the order of 10ms. In this case, the gains of the PI controllers are 
expressed in terms of machine parameters and response time as follows: 



(17) 



(18) 



(19) 



(20) 



(21) 



We used the method of compensation poles for its speed; it is clear that n 'not the only valid method for the 
synthesis of PI controller. 

MODEL OF THE MULTIPLIER MECHANICAL PART 

The mechanical party of the turbine includes three directional blades pitch and of length R. They are fixed to a 
drive shaft in a rotation speed H t , a multiplier Connected of gain G. This multiplier causes the electric generator. 
We can model all the three blades as one and the same mechanical system characterized by the sum of all the mechanical 
characteristics. Due to the aerodynamic blade design, we believe that the coefficient of friction with respect to the air is very 
small and can be neglected. Also, the speed of the turbine being very low, Friction losses will be negligible compared to the 
friction losses on the side of the generator. On the basis of these hypotheses, is then obtained a mechanical model consisting 
of two masses as shown in figure 6, the validity, relative to the complete model of the turbine, has already been 
verified(Lubosny. 2003; Usaola, Ledesma, Rodriguez, Fernadez, Beato, Iturbe, Wihelmi). 



Figure 6: Mechanical Model of the Wind Turbine 
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With: 

• J t : Moment of inertia of the turbine to the equivalent inertia of the three blades of the turbine. 

• J m : the moment of inertia of the DFIG. 

• / : The coefficient due to viscous friction of DFIG. 

• C m : mechanical torque on the shaft of the DFIG. 

• Jl m : the rotational speed of the DFIG. 

C em : the electromagnetic torque of the DFIG. 

By considering that the multiplier is ideal, That is to say, the mechanical losses are negligible, it is then modelled 
by the following two equations: 

C m = J (22) 

n m =G.n t (23) 

From Figure 6, we can write the fundamental equation of dynamics of the mechanical system of the mechanical 
shaft DFIG by: 

(^ + Jm)^ + f-n m = C m -C em (24) 
The total inertia J : 

y = (^+y m ) (25) 

CONTROL STRATEGIES OF THE TURBINE 

The control strategy is to adjust the torque appearing on the tree turbine speed so as to fix a reference. 
To achieve this, we will use a speed control (EL Aimani. 2004; Poitiers. 2003). 

According to equation (24) and (25): 

^ = j.(C n -f.l2 m -C m ) (26) 
The electromagnetic torque is: 

Cemjref PI • \Qref ^mec) (27) 
£l re f — G . f2 tur b i ne re f (28) 

The reference speed of the turbine corresponds to the optimum value corresponding to the speed ratio X pmax 
(fixed on the blade angle P to 2°). 

n - Xc Pmax v (?Q s 

iL turbine_ref R \^ y J 
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The couple thus determined by the controller is used as a reference torque of the turbine model as can be seen in 



figure 7. 
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Figure 7 : Block Diagram of the Controller with Mechanical Speed 

CORRECTION OF THE SPEED OF THE TURBINE 

Different technologies can be considered markers for control of the speed; in our case we describe two types of 
regulations: 

A. PI Controller with Anticipation 

We consider a correction proportional integral (PI): 

Cem_ref = (^1 + ~) ■ {pref ~ ^mec) (30) 

b-L : proportional gain and b 0 : integral gain, correction parameters are determined. It is necessary to increase the 
parameter b 0 to reduce the action of C„ wind torque. The natural frequency and damping coefficient are given by: 



f^o j er f+J+ b i "n 

Wn =J-and^ = — .- 



Can impose a time response and damping factor given was: 
b 0 = a)lJmdb 1 = 2 -^-f-J 

B. PI Controller with Phase Advance 

We consider a correction proportional integral (PI): 

a.i.s+a Q 



(31) 



(32) 



_ ai.s+a 0 / a 
u em_re/ T s+1 ' V re / mcc; 



(33) 



a 1; a 0 eti are the parameters to determine the corrector and s is the Laplace variable. It is necessary to increase the 
parameter a 0 to reduce the action of C„ wind torque. The natural frequency and damping coefficient are given by: 



gQ+/ , r _ bin 
J.T 311 ? - ag+f ■ 2 



(34) 
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SIMULATION RESULTS 



The simulation is done by imposing active and reactive power reference (P re f, Q re f), While the machine is entailed 
in variable speed 

Pref — V-Pm^pt' Qref ~ 0< 

with: 

r|: Performance DFIG; 

P m opt'- The optimal mechanical power 

The simulation results are given by the following figures. 
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Figure 8: Mechanical Speed with Anticipation Control 




Figure 9: Electrical Active Power Produced with Anticipation Control 
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Figure 10: Electrical Reactive Power Produced with Anticipation Control 



Index Copernicus Value: 3.0 - Articles can be sent to editor@impactjournals.us 



Correction of the Mechanical Speed for the DFIG Wind Turbine 



37 



phase advance 
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Figure 11: Mechanical Speed with Phase Advance Control 
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Figure 12: Electrical Active Power Produced with Phase Advance Control 
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Figure 13: Electrical Reactive Power Produced with Phase Advance Control 

CONCLUSIONS 

In our work, we have established the model of the machine using its power equations in the dq axis system related 
to synchronization. We have also developed the method of vector control power of the machine to know the order. 

The two methods of maximizing power were almost granny gives results (speed mechanical, power active and 
reactive) the only difference and the transient time in the case of PI with big anticipation and compared PI phase advance. 

APPENDICES 

• Nominal Power =1.5 (Mw) 
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• Stator Per Phase Resistance =0.0 1 2 (£2) 

• Rotor Per Phase Resistance=0.02 1 (£2) 

• Stator Leakage Inductance= 2.0372. 10" 004 (H) 

• Rotor Leakage Inductance= 1 .7507. 10" 004 (H) 

• Magnetizing Inductance= 0.0 1 35 (H) 

• Number Of Poles Pairs=2 

• Moment Of Inertia= 1 000 (Kg.M 2 ) 

• Friction Coefficient =0.0024 
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